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Abstract 

We study CP asymmetries in two-body decays of top squarks into neutrali- 
nos and sleptons at the LHC. These asymmetries are used to probe the CP 
phases possibly present in the stop and neutralino sector of the Minimal Super- 
symmetric Standard Model. Taking into account bounds from experimental 
electric dipole moment searches, we identify areas in the mSUGRA parameter 
space where CP asymmetries can be sizeable and discuss the feasibility of their 
observation at the LHC. As a result, potentially detectable CP asymmetries 
in stop decays at the LHC are found, motivating further detailed experimental 
studies for probing SUSY CP phases. 
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1 Introduction 



Supersymmetry (SUSY) [I] is a well motivated theory to extend the Standard 
Model (SM) of particle physics. SUSY models are not only favored by gauge cou- 
pling unification and naturalness considerations, but are also attractive from the 
cosmological point of view. For instance the lightest SUSY particle (LSP) is a good 
dark matter candidate if it is stable, massive and weakly interacting [2],[3]. Most in- 
terestingly SUSY models provide a number of new parameters, among some having 
physical phases which cause manifest CP violation [I] . In the SM, the CP phase in 
the quark mixing matrix, which is currently confirmed in B meson experiments 0E], 
cannot explain the observed baryon asymmetry of the universe [7j, and additional 
sources of CP violation in models beyond the SM are required [8]. 

In the Minimal Supersymmetric Standard Model (MSSM) [I], a set of remaining 
complex parameters is obtained, after absorbing non-physical phases by redefining 
particle fields. The complex parameters are conventionally chosen to be the Hig- 
gsino mass parameter the U(l) and SU(3) gaugino mass parameters M\ and M 3 , 
respectively, and the trilinear scalar coupling parameters Af of the third generation 
sfermions (/ = b, t, r), 

// = Me**", M 1 = \Mi\e i + 1 , M 3 = |M 3 |e^ 3 , A f = \A f \e*+ A f. (1) 

These phases contribute to the electric dipole moments (EDMs) of the electron, neu- 
tron and deuteron, which can be beyond their current experimental upper bounds [5J 
I9"HT2~|. However, the extent to which the EDMs can constrain the SUSY phases 
strongly depends on the considered model and its parameters [T314T8] . Thus mea- 
surements of SUSY CP observables outside the low energy EDM sector are necessary 
to independently determine or constrain the phases. 

For example the phases of the trilinear scalar coupling parameters A f have a sig- 
nificant impact on the MSSM Higgs sector [19]. Loop effects, dominantly mediated 
by third generation squarks, can generate large CP-violating scalar-pseudoscalar 
transitions among the neutral Higgs bosons j2U]. As a result the lightest Higgs bo- 
son with a mass of order 10 GeV or 45 GeV [2] cannot be excluded by measurements 
at LEP [22\. The fundamental properties and the phenomenology of CP- violating 
neutral Higgs boson mixings have been investigated in detail in the literature [23] . 

In addition, the phases of the coupling parameters Aj, )t)T in the third generation 
sfermion sector are rather unconstrained by the EDMs [T7l[T8]. Therefore it looks 
appealing to find the adequate CP observables in sfermion production and decays 
at high energy colliders like the LHC [21] or ILC [25], which are sensitive to the 
eventually large CP phases of A^ ttT . Third generation sfermions also have a rich 
phenomenology due to a sizable mixing of left and right states. 

The phases can drastically change SUSY particle masses, cross sections, branch- 
ing ratios [2"oT - |2"8"] . and longitudinal polarizations of final fermions [29]. Although 
such CP-even observables are very sensitive to the CP phases (the observables can 
change by an order of magnitude and more), CP-odd (T-odd) observables have to be 
measured for a direct evidence of CP violation. CP-odd observables are, for exam- 
ple, rate asymmetries of cross sections, distributions, and partial decay widths [30] . 
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However, these observables require the presence of absorptive phases, e.g. from 
loops, and thus usually do not exceed the size of 10%. 

Larger CP asymmetries in particle decay chains can be obtained with triple 
products of final particle momenta [31]. They already appear at tree level due to 
spin correlations. Triple product asymmetries have been intensively studied in the 
production and decay of neutralinos [3"2"H3T)| and charginos [3T)H3"8] at the ILC, also 
using transversely polarized beams [3H]- For recent reviews, see Ref. [10] . 

In this paper, we study CP asymmetries in two-body decays of a stop, 

tm^t + x-; m = l,2; i = 2,3,4; (2) 
followed by the subsequent two-body decay chains of the neutralino, 

X"^h+L 4^Xi + ^2; n = L,R, £ = e,fi, (3) 
and the top quark 

t^b + W- W^u e + £ 3 ; e = e,/i. (4) 

See Figure [1] for a schematic picture of the entire stop decay chain. Since the stop 
decay, Eq. (T5]), is a two-body decay of a scalar particle, only the spin-spin correlations 
of the outgoing fermions contain a CP-sensitive part. Spin-spin correlations are 
terms in the amplitude squared that include the polarization vectors of both the 
neutralino and the top. The angular distributions of their decay products are then 
correlated to each other due to total angular-momentum conservation. Thus for 
the measurement of the top-neutralino spin-spin correlations, it is vital to include 
momenta from both their decays into the triple products, otherwise the CP-sensitive 
information in that process will be lost [4"T] . 

Triple product asymmetries have been studied in two-body decays of stops [H] 
and sbottoms [12], however in their rest frame only. At colliders like the LHC the 
particles are highly boosted, which will generally reduce the triple product asym- 
metries [43JIII]. We also include the production of stops at the LHC, and analyze 
the stop boost distributions as well as their impact on the triple product asymme- 
tries in detail. In our numerical analysis we work within the minimal Supergravity 
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Table 1: Overview of sfermion two-body decays studied in the literature to analyze 
CP phases. Two-body decays are denoted by a — > be and three-body decays by 
a — > bed. The column 'CP Phases' lists the phases the process is sensitive to, and 
'Max ^4cp' gives the maximal asymmetries which have been found in the rest frame 
of the sfermion. Asymmetries in parentheses correspond to the decay W —> cs. 



Process 


Sub Decay 1 


Sub Decay 2 


CP Phases 


Max ^l C p 


Ref. 




t -> bW 

xi -> it 


W ->■ &/(cs) 
^X? 




40% 


en 


t-m ^ tXi 


t -> bW 

xl tlx\ 




(t>A t A^(t>M 1 


15% 


m 




t 

|_ 

xf -> ^ 

xf 

^ ^x? 


W ->■ &/(cs) 

^x? 

W ->■ fl/(cs) 


<t>A h ,4>!J, 


15(40)% 
5(10)% 
2(5)% 
2(5)% 


021 


f ffi 


x° -> ^x? 


Z ^fi 
Z — > cc 
Z ->■ 66 




5% 
15% 
20% 


133] 


q ->■ qx,i 


x° 




0Mi 


n/a 






x° 




0Mi 


15% 





(mSUGRA) framework of SUSY breaking [T], and explicitly add CP phases for the 
parameters A t , Mi, /x at the electroweak scale, fully taking into account the implica- 
tions on the EDMs. In addition, we calculate the amplitude squared for the entire 
stop decay in the spin-density matrix formalism [35]. The compact form of the am- 
plitude squared allows us to identify the optimal CP observables for squark decays. 
For a recent study of CP-odd observables in stop decays, including their production, 
but for neutralino three-body decays, see Ref. jl3]. A summary of the literature in 
sfermion two-body decays is shown in Table [TJ Three-body decays of stops have 
been studied in Ref. |4"o] . 

The paper is organized as follows. In Section [2J we identify the CP-sensitive 
parts in the amplitude squared, classify the CP asymmetries in top squark decays 
and discuss their dependence on the complex stop-top-neutralino couplings. Sec- 
tion |3] contains a discussion of the experimental EDM constraints on the SUSY CP 
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phases and asymmetries. For the LHC, we discuss in Section H] the SUSY parameter 
dependence of stop production, boost distribution and branching ratios, and give 
lower bounds on the required luminosities to observe the asymmetries over their 
statistical fluctuations. Section 2] also includes a discussion of angular and volume 
distributions of the triple products. We summarize and give our conclusions in 
Section In the Appendix, the stop and neutralino Lagrangians with complex cou- 
plings are reviewed, we give the phase space and calculate the stop decay amplitudes 
in the spin-density matrix formalism. We also show that absorptive phases can be 
eliminated from the T-odd asymmetries to obtain true CP asymmetries. 

2 CP Asymmetries in Top Squark Decays 

In this Section, we identify the CP-sensitive parts in the amplitude squared of the 
two-body decay of the top squark, see Eqs. ©-(jlj) and Figured! In order to probe 
these parts, we systematically classify T- and CP-odd asymmetries. These can be 
defined from various triple or epsilon products of the different final state particle 
three- or four-momenta Pt,Pb,Pe!,Ph,Pe 3 , respectively. Explicit expressions for the 
squared amplitude, Lagrangians, couplings, and phase space elements are summa- 
rized in the Appendix. 

2.1 T-odd Products 

The amplitude squared for the stop decay chain, see Figure [TJ can be decomposed 
into contributions from the top spin correlations, the neutralino spin correlations, 
the top-neutralino spin-spin correlations, and an unpolarized part, see Eq. (IE. 62}) . 
The amplitude squared is a sum of terms which are even and one term which is 
odd under a T transformation. Assuming CPT invariance, CP asymmetries are 
proportional to the T-odd term in the amplitude squared. Such a term originates 
only from the spin-spin correlations. It is given by the last summand in Eq. (1E.67j) . 
which is proportional to 

|T| 2 D \m{a t mi {b t mi )*}[pi,pt,P£ i ,P£^\, m = 1,2, 2 = 1,..., 4. (5) 

The left and right couplings a^ mi and b % mi are defined through the t m -t-j$ La- 
grangian |4Tj . 

=9* (<4 Pr + Pl) X° tm + h.c. , (6) 

which we give explicitly in Appendix O These couplings depend on the mixing in 
the stop and neutralino sector, and thus also on the CP phases <pA t , (p^ and (f>M x - 
The imaginary part of the coupling product, \m{a t rni {b t mj )*} , in Eq. ([5]) is multiplied 
by a T-odd epsilon product, for which we have used the short hand notation 

[pi , Pi , P3 , Pa] = pt P u 2 Ps p%, ( 7) 

with the convention £0123 = 1- Since each of the spatial components of the four- 
momenta changes sign under a naive time transformation, t — > —t, this product is 
T-odd. 
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In order to construct asymmetries that probe this CP-sensitive part, we consider 
an epsilon product of four four-momenta Pi 

£ = \Pi,P2,P3,Pa]- (8) 

The pi can in principle be any four linearly independent momenta in the two-body 
decay-chain of the top squark, Pi,Pt,Pb,Pe 1 ,Pe 2 iPe 3 i see Figured! Note however, that 
momenta from both the decay products of neutralino and the top have to be included 
in the epsilon product to obtain non- vanishing asymmetries [JTJ • Otherwise the top- 
neutralino spin-spin correlations are lost. Only they have CP-sensitive terms, since 
the stop decay, Eq. (j2J), is a two-body decay of a scalar particle^. 



2.2 T-odd Asymmetries 

For a given T-odd product £ , Eq. ([8]), we can now define the T-odd asymmetry of 
the partial stop decay width T |41j . 

= T(£ > 0) - T(£ < 0) = /Sign[g]|r| 2 dLips 

r(£>0) + r(£<0) f|T| 2 dLips ' [) 

with the amplitude squared |T| 2 , and the Lorentz invariant phase-space element 
(iLips, such that / |^| 2 rfLips = T. The T-odd asymmetry is also CP-odd, if 
absorptive phases (from higher order final-state interactions or finite-width effects) 
can be neglected [3T]. In Appendix El we also construct a CP asymmetry, where 
those contributions are eliminated. 

In general the largest asymmetries are obtained by using an epsilon product 
which matches the kinematic dependence of the CP-sensitive terms in the ampli- 
tude squared. In the literature, this technique is sometimes referred to optimal 
observables [Hj. Thus the largest asymmetry is obtained from the epsilon product 

£ = \phPt,pt 3 ,pti]- (io) 

Other combinations of momenta in general lead to smaller asymmetries, see Ref. |41j . 

Triple products of three spatial momenta can also be used to define asymme- 
tries [21] • In the stop rest frame, p~ = (m 4 -, 0), the epsilon product reads [UJ 

\PbPuPhiPh\ = m t Pt ■ (Pl 3 x P*J = rniT. (11) 

That triple product has been found to give the largest asymmetries in the stop rest 
frame, and other combinations of momenta for T lead to smaller asymmetries [41] . 
Note that the asymmetries of an epsilon product £ are by construction Lorentz 

3 If for example momenta from the top decay are not taken into account, only CP asymmetries 
can be obtained, which are sensitive to the CP phases (j)^ and 4>m 1 alone, which enter solely from 
the neutralino decay. Still in that case, a three-body neutralino decay is required [43l|44l|42] , or a 
two-body decay chain with an intermediate Z-boson |33) . However, the sfermion then merely serves 
as a production channel for neutralinos, whose CP properties are studied through their subsequent 
decays [331 H31 Ell S3- The asymmetries are then only of the order of 10%, which is typical for 
neutralino three-body decays , and also chargino three-body decays [35] , 
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Table 2: mSUGRA benchmark scenario, extended with SUSY CP phases at the 
electroweak scale. 



m 


m l/2 


A 


tan/3 4> At M (f) Ml 


70 GeV 


270 GeV 


500 GeV 


5 i?r 

5 



Table 3: SUSY parameters and masses at the electroweak scale for the benchmark 
scenario of Table [2] (£ = e, fi). 



Mi = 110 GeV 


M 2 = 206 GeV 


/i = 344 GeV 


\A t \ = 329 GeV 


m f = 472 GeV 


m- t2 = 612 GeV 


m~ bi = 560 GeV 


= 581 GeV 


m v o = 102 GeV 

X.i 


m v o = 189 GeV 


m v o = 350 GeV 


m x o = 376 GeV 


m + = 187 GeV 


m v ± = 376 GeV 

A2 


m- lR = 130 GeV 


m ih = 201 GeV 



invariant whereas those constructed with a triple product T are not [H] . The triple 
product asymmetries will therefore depend on the stop boost, = \pi\/E^, and are 
generally reduced if not evaluated in the stop rest frame |49]. We will discuss in 
detail the impact of the stop boost on the asymmetries at the LHC in Section HI 



2.3 Phase Dependence 

In order to analyze the phase dependence of the asymmetry A, Eq. OH]), we insert 
the explicit form of the amplitude squared in the spin-density formalism. As shown 
in Appendix [Fj we obtain 

A _ 2r] f Si g n (£)(Pb-Pu e ) [PhPt,Pi 3 ,P£i} dLips ^ 



m o — m- 



I(PfPe 3 )(Pb-Pu e ) dLips 



with the coupling function 



lm{aUbLY} 



(13) 



X 



Thus the asymmetry can be separated into a kinematical part, and the effective 
coupling factor 77, that is approximately independent of the particle masses, but 
governs the complete phase dependence of the stop-top-neutralino couplings. 

To discuss the phase dependence of rj (and A), we choose the mSUGRA frame- 
work of SUSY breaking pQ . We extend it by adding the CP phases 4>A t , and 4>m 1 
at the electroweak scale. We define a benchmark scenario in Table |2j and give the 
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Figure 2: Approximated CP-odd coupling factor \m{a\ 2 {b\ 2 )*} , Eq. (JHJ), (left), and 
complete effective coupling factor 77, Eq. (IT5|) . (right) for the decay t\ — >• tx°? as 
a function of the CP phases 4>A t an d </v All other phases are set to zero. The 
mSUGRA parameters are given in Table [2j 



other low energy parameters and masses in Table |3j The scenario is chosen to ap- 
proximately maximize the signals at the LHC, while still respecting current bounds 
from electric dipole moments, see the detailed discussion in Section [3J 

Within mSUGRA, the lightest neutralino is essentially a bino whereas the second 
lightest neutralino is a wino. For the stop decay via the second lightest neutralino, 
t\ — > tx®, the asymmetries will primarily depend on and M , through the influ- 
ence of stop mixing. We expect a weak dependence on (f> Ml , since x 2 nas a negligible 
bino component. In this limit of the neutralino mixing matrix [18], valid for a 
large part of the mSUGRA parameter space, the phase dependence of the CP-odd 
coupling factor ]m{a t mi (b t mi )*} for the decay t\ — > t%\ can be approximated as 



Im{e4(&i 2 )*} oc 



cos 2 9i 


m t \n 




tan f3 


f 


2 _ 


-Ml 



sin 



+ 



+ 



+ 



sin 2Bt 

sin20 t - 
tan/3 

sin 29i 

tan 2 /3 



m t M 2 



W - Mi 



m t 




M 2 


\fJL 


2 -Mi 



sm 



sin 



+ 



m t \n\ 



M 2 - Mi 



sm 



(<f> At + 20 M ) 



(14) 



with the stop mixing angle 9$ , defined in Eq. ( 1A.8|) . 

In Figure El we show the phase dependence of Im{a^ 2 (^ 2 )*} in its approximation, 
Eq. ( 1T4]) . for the mSUGRA scenario of Table [2j For comparison we also show the 
coupling factor rj, Eq. ( 1T3|) . without approximations. The deviations are due to 



S 



n/2 n 3n/2 2n n/2 n 3n/2 In 

4>A, <f>n 



Figure 3: Asymmetry Ait£\£z), Eq. (j9|, as a function of the CP phase 4>A t with 
(p^ = 0, Ti (left), and as a function of <^ with <pA t = 0, n (right), for the stop decay 
ti —> and X2 ~^ P-i&Ri see Figure [H in the stop rest frame. The mSUGRA 
parameters are given in Table [2J 

CP-even terms in the denominator of 77. The CP-even parts are generally functions 
of the cosine of the phases. Figure [2] illustrates that such terms have an important 
impact on the phase dependence of the CP asymmetries. 

We further illustrate this effect on the asymmetry A f|T2|) for the optimal epsilon 
product £ fTTU]) . for pp — >• tit\ production and the subsequent decays i x — )■ tyQ and 
X2 ~^ via the second lightest neutralino and a right slepton. In the stop rest 
frame, this asymmetry is equivalent to the triple product asymmetry with T = 
Pt ■ {pi 3 x Ph), as given in Eq. ( TTT|) . We use the short hand notation Aitiiis), and 
similarly for other asymmetries to indicate the momenta used for the triple product. 

In Figure El we show ^4(^1^3) as a function of <pA t (left panel), and (p^ (right 
panel). The dependence on <p Ml is numerically small as expected. We clearly see 
that the maxima of ^4(^1^3) ~ ±40% are not necessarily obtained for maximal CP 
phases (j>A t = tt/2, 3/27T. The reason is that the phase dependence of A(t£i£ 3 ) is 
almost governed by the coupling factor 77, see FigureEJ right panel. Since the CP-odd 
(CP-even) numerator (denominator) of rj has a sine-like (cosine-like) dependence on 
the phases, the maxima of 77, and thus in turn of A(t£i£z), are shifted away from <pA t , 
M = 7r/2,3/27r in FigureEJ Thus the asymmetry can be sizable for smaller values 
of the phases, which is favored by EDM constraints. For example, the asymmetry 
has a maximum at M w 3/87T for <pA t = 0, and for M = the asymmetry is 
maximal at <pA t ~ 1/57T. The positions of the maxima will remain when including 
the effects of the stop boost, as we will discuss in Section H~T1 This motivates the 
choice 4>A t = 1/57T, (p^ = at the weak scale in our benchmark scenario (Table [2]). 

3 Bounds from Electric Dipole Moments 

The CP-violating phases in SUSY models are constrained by limits on the T-odd 
electric dipole moments (EDMs) [T7j. The current experimental bounds on the 
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n/2 Ti 3n/2 2n 

<Pa, 



Figure 4: Electric dipole moment of Thallium as a function of 4>A t with <pA b = 
40/180 7T, <p Ar = 1/180 7T, 0m 3 = 10/180 tt and A/l = M = 0. The mSUGRA 
parameters are given in Table HJ 

EDMs of Thallium j9], neutron [10], and Mercury [TT] are, respectively, 

|d T1 | = 9 • 10~ 25 e cm, (15) 
|dn| = 3 ■ 10" 26 e cm, (16) 
Idjjgl = 3 • 10" 29 e cm. (17) 

A planned experiment for measuring the deuteron EDM aims at a sensitivity of [12] 

|d d | = 1 -3 • 10" 27 e cm. (18) 

These limits generically restrict the CP phases to be smaller than vr/10, in particular 
the phase (fi^ [15j[T8]. However, cancellations between different contributions to the 
EDMs can allow for larger CP phases of the order of 7r/4 [T7] . 

As an example of such cancellations, Figure |4] shows the EDM of Thallium as 
a function of <pA t , for (pA b = 40/180 ir, 4>a t = 1/180 ir, 0m 3 = 10/180 tx and 
0a/i = (fifi = 0. The other mSUGRA parameters are given in Table [2j At four 
particular values of <pA t , seen by distinct dips in the plot, cancellations reduce 
by more than two orders of magnitude. 

We now randomly scatter the phases 4>a t , 4>m 3 , 4>Mi an d 4>A t i setting the 
most constrained phase (fi^ = 0. For each scatter point, we calculate the maximal 
excess of the EDMs over the corresponding experimental limit with the software 
package CPsuperH [50]. For the set of scatter points, we then show the line of the 
minima of all the excesses obtained as a function of <pA t m Figure |5J For <pA t < ^ 
the minima are well below the experimental limit, and for <pA t > 71 om Y slightly 
above. So even large values of <pA t are not yet excluded by the current EDM limits. 
Note that a more exhaustive scan to find narrow cancellation gaps would probably 
make even smaller experimental EDM excesses possible. For comparison, we show 
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Figure 5: Experimental excess in electric dipole moments as a function of <f>A t for 
M = 0. The experimental excess is defined as the maximum max; \di((f>A t , •••)\/\df cp \, 
i = Hg, n, Tl, over all current experimental limits. The two curves correspond to 
the case where: all phases except <pA t are se t to zero ('</>a 4 only'), and the phases 
4>A h i 0a t , 4>Mi and (ft 

m 3 are randomly scattered ('all phases'). For a given value of 
(f>Ati the minimal excess is plotted. The mSUGRA parameters are given in Table [21 

in Figure [5] also the experimental excess in the case that all phases except <pA t are 
set to zero. This demonstrates that a conspiracy among the phases can considerably 
reduce the predicted electric dipole moments. 

We conclude that due to cancellations among different contributions to the EDMs 
only isolated points in the CP phase space can give large CP signals at the LHC. It is 
important to search for these signals, since the cancellations could be a consequence 
of a deeper model that correlates the phases. For instance, the existing EDM bounds 
can also be fulfilled by including lepton flavor violating couplings in the slepton 
sector [IB]. This is important when considering for example SUSY Seesaw models, 
where CP violation in the slepton sector is connected to the neutrino sector and 
Leptogenesis [5T1I52] . 

The interrelation between the EDMs and the SUSY CP phases show that CP 
observables have to be determined outside the low energy sector, for example by 
measurements at the LHC and ILC. Asymmetries of triple products generally depend 
in a different way on the SUSY parameters than the EDMs, in particular on the CP 
phases, giving rise to a potential synergy between experiments at the low and the 
high energy scale. Depending on the experimental results, the EDM bounds could 
be either verified, or the CP-violating sectors of the model have to be modified. 
CP-odd observables like triple product asymmetries would be the ideal tool for such 
independent measurements. 
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0.5 1.0 1.5 2.0 0.0 0.2 0.4 0.6 0.8 1.0 

m~ t [TeV] ft 



Figure 6: Total stop pair production cross section a(pp — > it*) at the LHC as a 
function of the stop mass (left). Normalized stop pair production distribution 
with respect to the stop boost factor /% (right). The leading order cross sections 
have been calculated at ^fs^p = 14 TeV using MadGraph [M]. 

4 CP Asymmetries at the LHC 

At the LHC, squarks are produced with a distinct boost distribution which reduces 
the triple product asymmetries in the laboratory (lab) frame. In contrast to the 
epsilon product asymmetries they are not Lorentz invariant. If the stop rest frame 
cannot be reconstructed, the triple product asymmetries have to be folded with the 
stop boost distribution. 

In this Section we will discuss the dependence of triple product asymmetries on 
the stop boost. To this end we first review the production cross section of stops pairs, 
show their boost distribution, and explain how the triple product asymmetries in the 
lab frame are obtained. After defining their theoretical statistical significance, we 
study the asymmetries, stop and neutralino branching ratios, and the significances 
in the mSUGRA parameter space. 

4.1 Stop Pair Production and Stop Boost Distribution 

The production of stop pairs at the LHC [53J 

p + p->t m + t* m , m = l,2, (19) 

dominantly proceeds via gluon fusion. As a result, the leading order cross section 
a{pp — > t m t*J is independent of any other SUSY model parameters than the stop 
mass [53]. The strong dependence can be seen in Figure [6] (left), where the cross 
section drops by six orders of magnitude with an increase of the stop mass from 
0.2 TeV to 2 TeV. For the LHC center-of-mass energy y/s^ = 14 TeV, we have 
calculated the hadronic cross section at leading order with the software package 
MadGraph [51] . Next to leading order (NLO) calculations for stop production have 
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0.0 0.2 0.4 0.6 0.8 1.0 n/2 n 3n/2 2n 

H fa, 

Figure 7: Asymmetries A{U x l z ), A{bhl z ), A{tbh), A{hl 2 h), see Eq. ©, for stop 
decay {t\ — > ix°, x° — >" ^-R, see Figure [T]) as a function of the stop boost (5^ (left). 
Asymmetry A{tl\l$) in the stop rest frame and the laboratory frame as a function 
of (f>A t with M = A/l = (right). The mSUGRA parameters are given in Table El 



been performed in Ref. [53], and soft gluon resummation for light-flavored squark- 
antisquark pairs at next to NLO has been performed in Ref. [55]. Mixed stop pairs 
can only be generated at order a A s with a strongly suppressed rate [53] . 
The produced stops have a distinct distribution in their boost 

Pi = ^ (20) 

along the direction of their momenta. In Figure M (right), we show the normalized 
boost distribution for several values for the stop mass. In our benchmark scenario 
(Table ED, the lightest stop has a mass of = 470 GeV and its' boost distribution 
peaks at (3^ ~ 0.9. Thus the decaying stop is highly boosted in the laboratory 
frame, which reduces the triple product asymmetries as shown in the next Section. 



4.2 Triple Product Asymmetries in the Laboratory Frame 

The asymmetries A, Eq. Q, which are based on epsilon products £ = [^1,^2,^3,^4], 
Eq. (|SD, are by construction Lorentz invariant. The corresponding asymmetries 
with triple products depend on the stop boost. In Figure [7] (left), we show the boost 
dependence of asymmetries for various triple product combinations T = (tiii^), 
(b£i£s), (tMi), and (^1^3) • The asymmetries are calculated for pp — > t\t\ produc- 
tion at the LHC with the subsequent decays t\ — > tyQ and X2 ~^ ^R- The functional 
dependence of the asymmetries on the stop boost is purely kinematic. In the stop 
rest frame, /3 f ~ = 0, they coincide with the corresponding epsilon product asymme- 
tries, and then decrease with increasing Note that the size of the asymmetries 
strongly depends on the choice of momenta, and largest values are obtained for the 
optimal triple product (^1^3) as given in Eq. 0. However this asymmetry also 
decreases most strongly with see Figure [7] (left). 
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A triple product asymmetry in the laboratory frame is then obtained by folding 
the boost dependent asymmetry A with the normalized stop boost distribution, 

A lah = - f % A(J3i) dft, (21) 
o Jo dPt 

with the production cross section cr = a{pp — > t m t* m ). This leads to a reduction to 
roughly 30-40% of the asymmetry A lah compared to A(fli = 0). This can be seen 
in Figure [7] (right), where we show both asymmetries as a function of the CP phase 
4>A t - Thus the stop boost effectively reduces the asymmetries, but does not change 
their shape with respect to the phase dependence. In the following we will focus on 
the folded asymmetry A(t£i£^) in the laboratory frame but similar discussions hold 
for asymmetries based on the other triple products. 



Theoretical Statistical Significance 

So far we have only discussed the asymmetry of the triple product, but in order to 
determine the statistical significance of a given asymmetry it is also necessary to 
take into account the total signal rate. For example, a large asymmetry does not 
necessarily result in a large significance if the signal rate is low as this does increase 
the statistical uncertainty. 

Assuming that the fluctuations of the signal rate are binomially distributed with 
the selection probability p = 1/2(^4+1), the significance of an asymmetry is 

S = - ^ v^Z, (22) 

with the integrated LHC luminosity C The cross section is 

a = F N x a(pp -»■ i m t* m ) x Br(t m -»■ tx°) x Br(t -> bW) x Br(W -»■ u e e) 

x F3r(x° ^ e + r n ) x Br(g~ xje"), (23) 

for m — 1,2, % — 2,3,4, and n = L,R. The combinatorial factor takes into 
account the possible W and neutralino x° decays into leptons with different flavors. 
We assume that the branching ratios do not depend on the flavor, i.e., Br(W^ — > 
v e e) = By(W — > i^/i), and Br(x° — > e + e~) = Br(x° — > fi + fl~). The factor is thus 
Fn = 8, if we sum the lepton flavors e, fi, and also the slepton charges, £~ and 

The statistical significance S is equal to the number of standard deviations that 
the asymmetry can be determined to be non-zero. For example a value of S = 1 
implies a measurement at the 68% confidence level. The minimal required luminosity 
is then 

c ^ H^ 1 )' (24) 

Our definition of the statistical significance S and the required luminosity L is 
purely based on the theoretical signal rate and its asymmetry. Detector efficiency 
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effects and contributions from CP-even backgrounds are neglected, which would 
reduce the effective asymmetry. The definitions have thus to be regarded as absolute 
upper bounds on the confidence levels and as absolute lower bounds on the minimal 
required luminosities, respectively. We will use these terms in the following to 
exclude mSUGRA parameter regions which cannot be probed at the LHC. In order 
to give realistic values of the statistical significances and required luminosities to 
observe a CP signal, a detailed experimental study is necessary. 

Also higher order corrections will be important and have to be included in a 
comprehensive analysis. Leading-order (LO) QCD calculations strongly depend on 
the a priori unknown renormalization and factorization scales. As a result, theoret- 
ical predictions of LO QCD cross sections and branching ratios are easily uncertain 
within a factor of 2, see for example Ref. [53]. However, we expect that asymmetries 
and distributions will have much smaller scale dependencies. Since the higher order 
corrections enter both in the numerator and the denominator, they might cancel par- 
tially, and may lead to a significantly reduced scale dependence of the asymmetries. 
For instance, the normalized stop boost distribution shown in Fig. |6] (right) only 
changes slightly if switching on initial and final state radiations in PYTHIA [56]. 
Soft gluon radiation could in principle change the spin structure of the stop decay 
process, see Figure [TJ However, only the top-neutralino spin-spin correlations are 
CP-sensitive, and thus they could only be altered by gluon radiation of the top 
alone. 

Finally, our asymmetries probe CP-violating phases of the stop-top-neutralino 
couplings, which only belong to the electroweak sector. However, we expect that also 
the effect of electroweak corrections to our observables and asymmetries is rather 
small. Although electroweak corrections to neutralino masses, for example, can be 
10% at one- loop level [57], and neutralino branching ratios for two-body decays may 
receive electroweak CP-even one-loop corrections of up to 16% in some cases [58] . 
we expect that asymmetries are less sensitive to these electroweak corrections, since 
they again enter both in the numerator and denominator. 

4.3 Phase Dependence 

We analyze the (f>A t an d <^ dependence of the triple product asymmetry A lah in the 
lab frame, as defined in Eq. f[2"Tj) . For stop pp — > iit\ production and subsequent 
decays ti — > ty^ and x\ ~~ >" ^Ri we show in Figure [8] (left) contour lines of A(t£i£ 3 ) 
in the 0^-0^ plane for our reference scenario (Table |2]). The maximal asymmetry 
in the lab frame is about 15%, whereas the dependence on (pM 1 is weak, due to 
the negligible bino admixture of \2 m mSUGRA scenarios. Note that the phase 
dependence of the asymmetry Figure M (left) closely resembles that of the effective 
coupling factor t], Eq. (IF. 79}) . shown in Figure [2] (right) for the same scenario. Thus 
the asymmetry directly probes the CP- violating structure of the stop-top-neutralino 
couplings. The kinematical part of the asymmetry is largely phase independent. 

To determine the potential observability of the asymmetry, we calculate the 
minimal required luminosity C, Eq. (|24|) . to observe a la deviation of the asymmetry 
above the statistical fluctuations. Contour lines of C in the 4>A t ~4>Li plane are shown 
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<f>A, <t>A, 

Figure 8: Triple product asymmetry A lab (t£i£ 3 ), see Eq. (J2Tp . in percent (left), and 
minimal bound on the luminosity £, Eq. (1211) . in fb _1 to detect such an asymmetry 
at lcr above fluctuations (right), as functions of the CP phases 4>A t , <fi^ w hh 4>m 1 = 0, 
for stop production at the LHC with J~s^, = 14 TeV and decay t\ — > i%2> X2 ~~ >" ^R 
see Figure [TJ The mSUGRA parameters are given in Table [2j 

in Figure (right). A large part of the parameter space can be probed with £ = 
500 fb" 1 . Even small phases close to the CP-conserving points <f>A t , 4>^ = 0, n can be 
probed with C = 50 fb" 1 . The significance also depends on the stop and neutralino 
branching ratios, see Eq. ( 123]) . which depend differently on the phases 4>A t , 4>^- Thus 
the contour lines of the luminosity do not directly follow those of the asymmetry. 

4.4 mSUGRA Parameter Space 

So far we have analyzed the asymmetries within our mSUGRA scenario (Table EJ). 
We now extend our study to the entire mSUGRA parameter space, and study the 
dependence of the asymmetries and the event rates on the parameters mo, m^, 
tan /?, and A Q . In such an mSUGRA framework, the production of light stops will 
dominate, pp — > tit*, followed by the subsequent decays ii — > tx% and X2 ~^ ^R- 

Dependence on tuq and m^j 

The asymmetries in the stop rest frame and in the LHC lab frame are rather inde- 
pendent of the masses of the stop and its decay products. Thus we find an almost 
constant asymmetry A lah {ti.il^) ~ 15% within the kinematical allowed region of the 
mo-mi/2 parameter plane. However the stop production rate and the branching 
ratios of the neutralinos sensitively depend on tuq and m^. 

In Figure [9] (left), we show contour lines of the leading order (LO) cross section 
a{pp — > tit\) x Br(ti — )■ tx°) i n the m -mi/ 2 plane. While the stop decay branch- 
ing ratio Br(ti — > ~ 0.1 is almost constant, the stop production cross section 
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Figure 9: Cross section a(pp — > tit{) x Br(ti — > tx®) (left), neutralino branching 
ratio into left slepton Br(%2 — > e+ &2) (middle), and X2 branching ratio into right 
slepton Br(x° — > e+ ^jj) (right), as a function of the mSUGRA parameters m and 
mi/2- The other mSUGRA parameters and the CP phases are given in Table [21 The 
dark red area is excluded by direct SUSY searches. 

strongly decreases with increasing stop mass, see Figure M (left). Therefore the cross 
section also rapidly drops with increasing m^, which clearly can be seen in Fig- 
ure [9J (left). The dashed lines m x o = m e - Hi indicate the neutralino decay thresholds 
into sleptons. We show their branching ratios Br(%2 — > e+ &i) anc l B r (X2 — > e+ &li) 
in Figure M (middle) and (right), respectively. The branching ratio into right se- 
lectrons reaches up to 12% for small values of mo and mi/2- The decay fraction is 
considerably reduced for mi/ 2 > 300 GeV, where decays into the Z and the lightest 
neutral Higgs boson open, see the dashed line m x o — m x o = rriz in Figure [9j The 
significance for observing an asymmetry for the decay X2 ~^ is suppressed to 
that for the decay X2 ~^ f° r most of the mo-mi/2 parameter space. This is 
because the branching ratio Br(x° — > ^l) is larger for higher values of m^, where 
however the stop production cross section is suppressed. We thus find that both m 
and m.1/2 should be as small as possible to obtain the highest process rates. 

These findings are confirmed when we calculate the minimal luminosity C, Eq. (]24[) . 
which is required to observe the asymmetry at la for stop production pp — > tit^ at 
the LHC, with the subsequent decays t\ — > anc l X2 ~^ ^l(r), see Figure [TUl (left) 
and (right), respectively. The mSUGRA parameters and the CP phases are given 
in Table [2j For right sleptons, the required integrated luminosity can be as low as 
2 fb _1 for m ~ 50 — 100 GeV and m^j ~ 200 GeV. Asymmetries involving left 
sleptons require at least C ~ 150 fbr 1 for m rs 50 GeV and mi/ 2 ~ 340 — 420 GeV. 
In the overlap between the allowed kinematical regions for decays into left and right 
sleptons, the rates can be comparable, though. As the two asymmetries have the 
same magnitude but opposite signs, A eL (t£i£ 3 ) = —A iR {tiii 3 ), this would lead to a 
reduction of the total asymmetry in this intermediate region. In Figure [TUl we also 
indicate the mQ-m-1/2 parameter region where three-body decays X2 ~^ ^ + ^~Xi have 
been used to study CP asymmetries [4"3] . 
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Figure 10: Contour lines in the rriQ-mi/2 plane of the minimal required luminosity 
£, Eq. (1231) . to observe the asymmetry A lah {t£i£ 3 ), Eq. (]2"Tj) . at la above fluctuations 
for stop production at the LHC with y/s^ = 14 TeV and decay t± —> tx®, via a left 

slepton X2 ~~ >" ti&L (left) and right slepton x\ ~ * (right), see Figured] The 
other mSUGRA parameters and the CP phases are given in Table El The dark red 
area is excluded by direct SUSY searches. 

Dependence on Aq and tan/3 

In Figure mj we show contour lines in the Ao-tan /3 plane of the asymmetry A{ti\^) 
for stop production, pp — > t\t\, (and decay t\ —> tx®, X2 ~~ >" ^R-> see Figured]) and 
the corresponding minimal luminosity to observe it at la above fluctuations. All 
other parameters are given as in the reference scenario, Table HJ The asymmetry is 
essentially independent of Aq and tan/3 over the region of interest, slightly peaking 
at 16.5% at Aq ~ 500 GeV and tan/3 « 10 — 15. For small values of tan/3 < 2, a 
reduced stop mixing leads to smaller asymmetries. On the other hand, an increasing 
value of A Q also enhances the stau neutralino branching ratios, X2 ~^ ^ s a result 
of this, the dependence of the required luminosity on A and tan (3 is dominated by 
the rate, and the lowest luminosities are achieved for small values of Aq and tan/3, 
see Figure dH (right). 

4.5 Angular and Volume Distributions of Triple Products 

A measurement of an asymmetry requires the reconstruction of the spatial momenta 
of the three particles which form the triple product. It has to be determined whether 
the three momenta form a left- or right-handed trihedron. This can certainly only be 
answered by a detailed experimental simulation for the LHC. However to gain some 
insight, we briefly discuss the distributions of the stop decay width with respect to 
the volume and opening angles of the trihedron spanned by the momenta. 
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Figure 11: Contour lines in the A -tan/3 plane of the asymmetry ^4. lat> (t£ 1 £ 3 ) , 
Eq. (1211) . in the laboratory frame (left), and the minimal required luminosity C, 
Eq. (j24j) . to observe it at la above fluctuations (right), for stop production at the 
LHC with y/s^ = 14 TeV and decay i\ — > tx®, X2 ~~ >* h^R, see Figured] The other 
mSUGRA parameters and CP phases are given in Table [2j The dark red area is 
excluded by direct SUSY searches. 
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Figure 12: Normalized distributions of the stop decay width T with respect to the 
cosine of the three opening angles fa {i = t, £1, £ 3 ) of the parallelepiped (p t , p^, Pe 3 ), 
e.g. fa = < (pt, Pi! x Pe 3 )- Results are shown for stop boost factors 0^ = (left), and 
fit = 0.8 (right). The mSUGRA parameters and CP phases are given in Table [2J 

In Figure [T2J we show the distribution of the stop decay width, dT /d cos (pi, 
with respect to one of the three opening angles of the trihedron formed by the 
momenta p t , p^, Pi 3 , where fa (2 = t, £3) is the opening angle between one particle 
momentum and the plane formed by the other two, e.g., fa = <(pt, pe 1 x p^ 3 ). In 
Figure IT21 (left), the distribution is shown in the stop rest frame, and in Figure IT21 
(right) for a boost factor (3^ = 0.8. The parallelepiped collapses to a plane for 
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Figure 13: Normalized distributions of the stop decay width T with respect to the 
normalized volume V = pt ■ (pe 1 x P4)/ (|pt| |p^ | IP4 1) of the three particle momenta 
which form the parallelepiped. Results are shown for stop boost factors (3^ = and 
fit = 0.8. The mSUGRA parameters and CP phases are given in Tabled 

cos (f>i = 0. For angles close to this configuration it will be difficult to reconstruct 
the handedness of the trihedron. For larger stop boost factors, the distributions 
become even more symmetric and the peaks move closer to cos<fti = 0. Only for 
4>h = ^(Pii)P^3 x P*)' Figure [T21 (right), double peaks are formed at « ±45°, 
which should simplify the determination of the handedness in a boosted frame. 

In Figure fT3| we show the normalized distribution dT/dV with respect to the 
normalized volume V = £> t ■ (p^ x p£ 3 ) of the unit vectors p = p/|p| for a stop 
boost factor f3^ = and (3^ = 0.8. This succinctly demonstrates that the event 
configurations mainly have the form of a collapsed parallelepiped, in particular for 
large stop boosts. 

5 Summary and Conclusions 

We have analyzed observables of CP violation in the two-body decays of a light stop 

ti -)• t + xl 

The CP-sensitive parts appear only in the top neutralino spin-spin correlations, 
which can be probed by the subsequent decays 

t ->■ b + W, W^u e + £ 3 

This specific decay chain of the light stop was found to give the largest signal of CP 
violation over most of the mSUGRA parameter space. 

Due to angular momentum conservation, the decay distributions of the final 
state momenta are correlated to each other. Asymmetries of triple products of 
three spatial momenta, as well as epsilon products of four space-time momenta are 
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ideal CP observables. They are sensitive to the CP phases of the trilinear coupling 
parameter A t , the higgsino mass parameter /i, and the gaugino mass parameter Mi. 

We have analyzed the asymmetries and event rates in an mSUGRA framework, 
with CP-violating phases explicitly added at the weak scale. For a strong stop 
mixing, the asymmetries are rather constant in the mSUGRA space, and reach up 
to 40% in the rest frame of the stop for 0^ = 1/57T and ^ = 0. The influence of 
<f>Mi is negligible for a wino-like neutralino. 

The asymmetries of triple products are not Lorentz invariant and therefore frame- 
dependent. Since the stops are produced highly boosted at the LHC, the asymme- 
tries are reduced by a factor of about three when evaluated in the laboratory frame, 
compared to the stop rest frame. Luminosities of at least 10 fb~ for small values of 
m i/2 %j 300 GeV, m < 150 GeV and tan/3 ps 5 are required to observe a CP signal 
of la above statistical fluctuations at the LHC. 

There are already stringent constraints on SUSY CP phases from electric dipole 
moments, but large phases and CP asymmetries at the LHC are possible when taking 
into account cancellations between EDM contributions from different phases. Since 
the asymmetries of T-odd products depend differently than EDMs on the SUSY 
parameters, in particular on the CP phases, asymmetries would be an ideal tool for 
independent measurements at the LHC and ILC. 

Clearly, the measurability of the asymmetries and thus of the CP phases can 
only be addressed properly in a detailed experimental analysis, which should take 
into account background processes, detector simulations and event reconstruction 
efficiencies. We want to underline the need for such a thorough analysis, to explore 
the potential of the LHC to probe SUSY CP violation. 
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Appendix 
A Stop Mixing 

The masses and couplings of the stops follow from their mass matrix 

m| e~ l ^ l m t \Ki\ 
?m = ~(tl, t* R ) I I , (A.l) 



e^m t |A^| m- 



2 

t R 



with 



m\ L = M? + Q- jjsin 2 ^m|cos2/3 + m 2 , (A.2) 
2 

m~ R = M~ + ~m| sin 2 W cos 2/3 + m 2 , (A.3) 

with the soft SUSY-breaking parameters Mq, My, the ratio tan/3 = ^ of the 
vacuum expectation values of the two neutral Higgs fields, the weak mixing angle 
9 W , the mass mz of the Z boson, and the mass m t of the top quark. The CP phase 
of the stop sector is 

<k = axg[Af], (A.4) 

A t - = A t -fi*catp, (A. 5) 

with the complex trilinear scalar coupling parameter A t , and the higgsino mass 
parameter /x. The stop mass eigenstates 

*1 A ~,t ( l L 



w v &■ 1 ' (A ' 6) 



are given by the diagonalization matrix 



e i<f>i cos s j n q. 

ft* = I I , (a.7) 



sin 9f e ^* cos 9i 



with the stop mixing angle 9$ 

-m t \M 



1 2 

mj — m 



cos% = "" rr " sin% = , tL = (A. 8) 

y m ? I A t I 2 + (™| - m t ? L ) 2 yj rn 2 t\ A i\ 2 + 0| - m 2 L ) 2 



The mass eigenvalues are 

2 1 

m,~ = - 

tl,2 2 



m~ + m~)T\/(ra 2 - - m| ) 2 + 4m 2 |A t -| 2 . (A.9) 



Note that for |At| S> cot /3 we have 
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B Neutralino Mixing 

The complex symmetric mass matrix of the neutralinos in the photino, zino, higgsino 
basis (7, Z, H®), is given by 



( M 2 sl + M lC l {M 2 -M 1 )s w c w \ 

(M 2 - MO M 2 c 2 w + Mi 4 m z 

m z iis 2 p -/ic 2/3 



[B.10) 



with the short hand notation for the angles s w = sin^, c w = cos9 w , and s 2 p = 
sin(2/3), c 2/ 3 = cos(2/3), and the SU{2) gaugino mass parameter M 2 . The phases of 
the complex parameters Mi = |Mi|e"^ M i and fi = |/i|e l<?ifl can lead to CP- violating 
effects in the neutralino system. We diagonalize the neutralino mass matrix with a 
complex, unitary 4x4 matrix N [I], 

N* ■ M x o ■ N 1 = diag(m x o, . . . ,m x o), (B.ll) 

with the real neutralino masses < m^o < m^o < m^o < m^o . This can be achieved 

Xl A2 A3 X4 

by diagonalizing the hermitian matrix -M^ ■ -M x o with a unitary matrix U 

U ■ {M[ ■ M x o) ■ C/t = diag(m 2 ? , . . . ,m 2 x o). (B.12) 
The neutralino mixing matrix N can then be expressed as 

N = diag(e^ l/2 , . . . , e^ 4/2 ) • U, (B.13) 
with the phases fa of the elements of the diagonal matrix U* ■ Ai x o ■ U\ 

fa = arg(ET • M x o ■ U%, % = 1, . . . , 4. (B.14) 

C Lagrangians and Couplings 

The interaction Lagrangian for stop decay t m — > tXi is given by [4"T] 



^0 = 5 t (a^ P fl + 6^ P L ) x° i~m + h.c, (C.15) 

with Pl,r = (1 T 7s)/2, an d the weak coupling constant g = e/ sin6 w , e > 0. The 
couplings are 

2 2 

a mi = ^ ] i'R-mn) ^m> ^mi = ^ ] (^mJ ^in' (C.16) 

n=l n=l 

with the stop diagonalization matrix 7£*, see Eq. (I A. 7ft . and 
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with, in the photino, zino, higgsino basis (7, Z, H®, 



f L 

J tt 



-y/2 



1 (\ 2 2 \ 2 

— ---sin 6 W ) N i2 + - sin 6 w N a 

cos U w \ 2 6 I 6 



R 



■hi 



h L 

lb ti 



Y, 



2 sm9 w (ttm0 w N* 2 - N^] 



3 

(h*y = Y t (N* 3 smf3-N* 4 cos/3), 



v/2 mw sin /3 

and the mass of the W boson. 

The interaction Lagrangian for neutralino decay x° — > followed by 

X°£ ± with I — e, 11 is given by 



(C.18) 

(C.19) 
(C.20) 
(C.21) 



R,L 



with the couplings 

fli = V2 



sin 2 9 W N i2 + sin 9 w N a 



cos6^ \2 
3l = V2sme w {taiiO w N; 2 -iV*). 
The Lagrangian for I = r with complex A T are given, for example, in Ref 



(C.22) 

(C.23) 
(C.24) 



D Kinematics and Phase Space 

For the squark decay t m — > tx®, we choose a coordinate frame in the laboratory 
(lab) frame such that the momentum of the top squark t m points in the ^-direction 



LI 

Pt 



(E h 0, 0, |p f |), 

(E t , \p t \ sin6 t , 0, |p t | cos6> t ). 



The decay angle 6 t = <X(pt, Pt) of the top quark is constrained by sin #™ ax = 
for |pf| > |p~| = A^(m?,m 2 ,m 2 )/2mf, with the triangle function A(a, b, c) 
b 2 + c 2 — 2(ab + ac + be). In this case there are two solutions 



(D.25) 
(D.26) 

ip;-i/iPfi 



ipfi 



(m~ + m 2 — m 2 )|p t -| cos 6 t ± E{</\(m 2 , m 2 , ni^o) — MPt\ 2 m 2 sin 2 6 t 



2 1 p t 1 2 sin 2 0* + 2m 2 



(D.27) 



For |pf| < |p~| the angle 9 t is unbounded, and only the physical solution |p+| is left. 
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The momenta of the subsequent decays of the neutralino Xi ~ > h^] £ Xi^2 (ED? 
and those of the top quark t — > bW, W — > £ 3 U£ (jlj), can be parametrized by 

p% = E b (l, sin 9b cos 4>b, sin 9b sin (p b , cos 9b), (D.28) 

p% = ^(1, sin cos 0i, sin sin 0i, cos ^i), (D.29) 

Ee 2 (l, sin ^2 cos 02, sin 6 , 2 sin0 2 , cos #2), (D.30) 

p% 3 = ^(1, sin 6» 3 cos 03, sin ^3 sin 03, cos 6» 3 ), (D.31) 

with the energies [36] 

m 2 — mj mj — m 2 

£l ~ 2(E x o - |p x o| cos9 Dl ) ' 12 ~ 2(E i - cos9 D2 ) ' 1 ' ] 

2 2 2 

771 _ m t ~ m Vi/ 771 _ m W (y) oq\ 

£ ' 6 ~~ o7^ i i a — V' e 'i ~ o7^ i i a — V' \u.66) 

2(E t - \pt\ cos9 Db ) 2(E W - \p w \ cos9 D3 ) 

and the decay angles, 9 Dl = <(p x o, p £l ), 9 D2 = <{pi,Pe 2 ), 9 Dz = <(pw,Pe 3 ), and 
®D h = <(pt,Pb), for example, 

cos9 D -. = : : — , (D.34) 

IP* — p*l 

with the unit momentum vector p = p/|p|. The other momenta p x o, p^, and pw, 
follow from momentum conservation. 

The Lorentz invariant phase-space element for the squark decay chain, see Eqs. fl2])- 
(J4]), can be decomposed into two-body phase-space elements [36|l60] 

dLips(s { ;p n ,p x o,p b ,p h ,p h ,pe 3 ) = -^-^ ^ dLips(s f ;p t ,p x o) 

xds x o dLips(s x o;p£ 1 ,pi) ds| dLips(sf,p x o,p£ 2 ) 

xds t dLips(s t ;pb,Pw) ds w dLips(s w ;pe 3 ,p Vi ), (D.35) 

where we have to sum the two solutions \pf | of the top quark momentum, see 
Eq. (1D.27j) . if the decay angle 9t is constrained. The different factors are 

I ip 12 

dLvos(si\p t ,p x o) = — | - i- | _ \ 1 1 sin0 t d9 t , (D.36) 

dLips(s x o;pi 1 ,pi) = n/ ^ 2 ^Ji^ ^,2 (D.37) 



1 


|P*| 2 




8tt \E t 


Pf cos9 t - 


-E~ t \p t \\ 


1 


Ip^iI 2 


dQi, 


2(2tt) 2 


m x o - mj 


1 


IPfel 2 


dQ 2 , 


2(2tt) 2 


mj - m x? 


1 


In I 2 
IP&I 


dflb, 


2(2vr) 2 


m\ - m 2 w 



dLips(sfp x o,pi 2 ) = nM ^ 2 — 2 — 2 —^d£li-, (D.38) 



dLips(s t ; Pb, Pw) = , N2 — , o - (D.39) 



I P«J 



2 



dLips{s w ]p h ,p Vt ) = — dn 3 , (D.40) 
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with Sj = p 2 , and dfL = sin#~ d9j d<f>j. We use the narrow width approximation 



m\ 2 ds 3 = (D.41) 
riijL j 



for the propagators 



Sj — + irrijL j 



which is justified for Vj/rrij <C 1, which holds in our case for particle widths 
Tj < C(lGeV), and masses rrij ~ (D(lOOGeV). Note, however, that the naive 
0(T/m) -expect at ion of the error can easily receive large off-shell corrections of an 
order of magnitude and more, in particular at threshold, or due to interferences with 
other resonant or non-resonant processes. For a recent discussion of these issues, 
see, for example, Ref. [61]. 



E Density Matrix Formalism 

The amplitude squared for the entire stop decay chain, Eqs. ([2])-(j4]), has been calcu- 
lated in Ref. [41], by using the spin formalism of Kawasaki, Shirafuji and Tsai |62j . 
We calculate the amplitude squared in the spin-density matrix formalism [4"5|l63]. 
which allows a separation of the amplitude squared into contributions from spin 
correlations, spin-spin correlations and from the unpolarized part. In that way, the 
CP-sensitive parts of the amplitude squared can be easily separated and identi- 
fied. The stop decay amplitude squared in this formalism, however for a subsequent 
three-body decay of the neutralino, was calculated in Ref. [43] . 

In the spin-density matrix formalism of Ref. [IS], the amplitude squared of the 
stop decay chain, Eqs. (JS])-®, can be written as 

|T| 2 =|A(t)| 2 |A0a| 2 |A(W0| 2 |A(£)| 2 x 
E P*>$)%K PnM) KXi D 2 (£) PDi {t) x $ PD4 (W)x kK . (E.43) 

The amplitude squared is composed of the propagators A(j), Eq. (1D.42[) . of particle 
3 = t] Xi\ W or £, and the un-normalized spin density matrices Adi(X;), 
PD 3 (t), and Pd a (W), with the helicity indices Aj, A^ of the neutralino, the helicity 
indices X t , X' t of the top quark, and those of the W boson, A&, X' k . The decay factor 
for the selectron is denoted by D 2 (t). 

The density matrices can be expanded in terms of the Pauli matrices 

+<a; < K (E-44) 
PdM)^ = Kx>Di + <T b XiK Z b Dl , (E.45) 
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h>x t + a a x , Xt 

XI* 



with an implicit sum over a,b = 1,2,3. 

The polarization vectors of the W boson fulfil the completeness relation 

PW,fi Pw,v 



(E.46) 
(E.47) 



E 



-9v» + 



(E.48) 



ir 



with Pw £ ^ k — 0- Similarly the spin four-vectors sf, a = 1,2, 3, for the top quark t, 
and s fe , b — 1, 2, 3, for the neutralino also fulfil completeness relations 



p?p u t 



m 



2 ' 



EJ>,tM b,U 
S x° S x° 



p^p% 

2 ' 
Xi 



(E.49) 



and they form an orthonormal set 

s a t -s b t = -5 a \ s a t -p t = 0, s a x o-s b x o = -5 ab , s a x o-p x o = 0, (E.50) 

with the notation p^ = p^/m. 

The expansion coefficients of the matrices (1E.44I) - (1E.46I) are 



D 



D 



snab 
D 



9 2 
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(\a 




9 2 


(+) 


2 




9 2 


(+) 
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9 2 
2 


(I. 



t 1 2 i \it 1 2 



\ h mi\ ) (Pt ■ P x °) ~ 9 Re{a mi (b mi )*}m t m x o 



A, ' 



J 1 2 \ut |2 



l^r)^(Pxr<) 



J |2 lit |2 



l 6 mil ) m x?(^" S x°) 



(E.51) 

(E.52) 
(E.53) 



t |2 



I6 1 



t \2 \ ( a „b 



si ■ s x o)m t m x o 



+g 2 Re{al l (bl l )*} (s a t ■ p x o)(4 9 ■ Pt) - « ■ s*o)(p x o • p t ) 



^3 



-g 2 iMaLi(Kni)*}[4, pt, s b x o, Px o] 

^l/^| 2 Ko-mJ), 

(V 2 l/#r^(4r^)' 



(E.54) 

(E.55) 
(E.56) 

(E.57) 



(+) 



ym* + " (P& " OfH (l^^a Pb,p}> (E.58) 



(E.59) 
(E.60) 
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with the couplings as defined in Appendix £3 and the short hand notation 

\P1,P2,P3,P4\ = £»vap Pi P2 Ps £\ ^0123 = 1- (E.61) 

The coefficients for squark decay, Eqs. flE.5ip - flE.54p . agree with those given in 
Ref. [12] • The signs in parentheses hold for the charge conjugated processes, that 
is t* m ->■ tx Q i in Eqs. flE32j) and (JESHJ), x°i ~> in Eq. flETM]) . t -> 6W - in 

Eqs. ( ESZj) and E5g) . and finally W~ -> z/^J in Eq. fiE59|) . 

Inserting the density matrices, Eqs. f IE . 44p - (IE . 47 j) into Eq. (IE.43j) . we obtain 

|T| 2 = 4|A(t)| 2 |A(x°)| 2 |A(iy)| 2 |A(£)| 2 x 

[D D x D p / + D x Y> a D £°f + S 6 Di D p / + E b Dl E^f ] /J 4pCT D 2 , (E.62) 

with an implicit sum over a, b. The amplitude squared \T\ 2 is now composed into an 
unpolarized part (first summand), into the spin correlations of the top (second sum- 
mand), those of the neutralino (third summand), and into the spin-spin correlations 
of top and neutralino (fourth summand), in Eq. (IE. 621) . 

With the completeness relation for the W polarization vectors (IE.48j) . we find 



D^ 7 Di pa 

spa pa n 
^DZ ^4 per 



2g 4 (PfPe i )(Pb-Pv e ), 
( +)2m^ 4 (sj • pt 3 )(p b ■ Pu t ), 



(E.63) 
(E.64) 



and the sign in parenthesis for the charge conjugated decay, t — > bW~, W~ — > vit~^. 
By also using the completeness relations for the top and neutralino spin vectors, 
Eq. (IE.49p . the products in Eq. (IE.62[) can be written as 



spa spa pa j~, 

^4 pa 



9 6 {\aL\ 2 



(Pb-Pu e ) 



(P x 9 ■ Pt){Pi 3 ■ Pt) - m 2 t (Px? ■ Pis) 



x, 



(E.65) 



spb spb 



spab spb spa pa t-\ 

^Di ^Ds u 4pa 



C-)2 



— a 



t 1 2 \it 1 2 



i-R\2 



X 



m 2 Y o(pt -Ph)- (p x ° ■ Pt)(Ph ■ P x o) 



(E.66) 



+ (Pt ■ Ph){Pt ■ Ph) m \] ~ (Pti ■ Ph) m l° m t 
- X ° ■ Pt ) (Pt ■ Pta ) (P x o - Pti) (Pb- Pu t ) 
2g 8 Re{ai i (bi i y}\f*\ 2 m x om t (p b ■ p Ut ) 

(P x ° ■ Pt)(pi! -Ptz)- (p x ° ■ Pe 3 )(Pt ■ Ph) 



x 



{ t ) 2g 8 lm{al ni (b t mi )*}\f^\ 2 m x om t ( Pb ■ p vl )\p^PuPh , P x f\ ■ 

(E.67) 
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Note that \pi 31 PtiPhiP x ] = \phPt,Pi 3 ,P£j] due to momentum conservation, with the 
short hand notation Eq. (IE.61[) . The signs in parentheses in Eqs. (1E.66[) and f]E.67|) 
hold for the charge conjugated process i — > bW~, W~ — > v^. In order to obtain 
the terms for the decay into a positively charged slepton, x° — \ followed by 

— > I2X11 one nas to reverse the signs of Eqs. (1E.66I) and (1E.67|) . This is due to 
the sign change of the term E^, , see Eq. (IE.56I) . 



For the decay into a left slepton \°i ~ > ^i ^li ^qs. AE.55j) . (1E.56j) . and (1E.60j) read 



Di = y|^| 2 Ko-m|), (E.68) 



S 6 n. = ^gVlfm^-pe,), (E.69) 



,2 



D^) = 9Vkr(rnj-m' x o), (E.70) 

respectively. The expressions for Eqs. (IE. 661) and (1E.67|) have to be changed accord- 
ingly. The sign in parenthesis in Eq. (IE. 691) holds for the charge conjugated process 



F Stop Decay Widths and Asymmetry 

The partial decay width for the stop decay t m — > ty® in the stop rest frame is 



ntm -+ 1%) = y- — 4nml *- d, (F.71) 

with the decay function D given in Eqs. (IE.51I) . For the decay t m —> b"Xj we have [26] 

-+ bxt) = mnm ^ g 2 (ifc^l 2 + 14/) , (F.72) 

whith the approximation m;, = 0. The stop-bottom-chargino couplings are [26] 

4' = ^ti^n, 4' = + K2 V j2 Y u (F.73) 



with the stop diagonalization matrix TZ l , Eq. (1A.7j) . the diagonalization matrices U, 
V for the chargino matrix pQ, 

U* ■ M x ± ■ y f = diag(m x ±, m x ±), (F.74) 
and the Yukawa couplings Y t , Eq. ( 1C21 1) . and 



n = To ^ ■ (F.75) 
V2 mw cosp 

The stop decay width for the entire decay chain, Eqs. (EJ)- (J4j) , is given by 
r(t ->■ vixlblihh) = y l T | 2rfLi P s ( s < ',Pu e ,P x o,Pb,Pi^Pi2,Pi 3 ), (F-76) 
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with the phase-space element c/Lips given in Eq. f]D.35|) . 

We obtain an explicit expression for the asymmetry, if we insert the amplitude 
squared |T| 2 flR62j) into Eq. ©, 

A = I Sign(g) gT d *p° ^2 ciLips 

/ D D x D 4pr7 D 2 dHps ' [ ' 

where we have already used the narrow width approximation of the propagators, see 
Eq. (ID. 41}) . In the numerator, only the spin-spin terms of the amplitude squared 



remain, since only they contain the epsilon product S, see Eq. (JSJ). The other terms 
vanish due to the phase space integration over the sign of the epsilon product, 
Sign(£). In the denominator, all spin and spin-spin correlation terms vanish, and 
only the spin-independent parts contribute. Inserting now the explicit expressions 
of the terms Eqs. (1E.51I) - (1E.67|) into the formula for the asymmetry, Eq. (1F.77j) . we 
find 

A o /Sign(5)(p 6 -^J \pi,Pt,Ph>Pix\ dU P s ,p 7 m 

A = 2 " K„-ml)/(p ( . P „)(p t . P „ ( )< i Lip S ' (R78) 



with the coupling function 



Im{<i( & mi)*} 



\ (l«Ll 2 + I^Ll 2 ) 2m x ; mt - M«L(&L)*> 

Xi 



(F.79) 



G Contributions from Absorptive Phases 

The T-odd asymmetries A, as defined in Eq. (Q, are not only sensitive to the CP 
phases in the stop-neutralino system. Absorptive phases could also enter from final- 
state interactions or finite- width effects. Although such absorptive contributions are 
a higher order effect, and thus expected to be small, they do not signal CP violation. 
However, they can be eliminated in the CP asymmetry 

A CP = ^(A-A), (G.80) 
where A is the asymmetry for the stop decay, 

i m -H + Xil with t^b + W + ; W + -H>i + 4; (G.81) 

+ + (G.82) 

and A the corresponding asymmetry for the decay of the charge conjugated stop, 

C-^+Xi! with t^b + W~; W-^vt + l-, (G.83) 

and the same neutralino x° decay as in Eq. (1G-.82j) . In forming the CP asymmetry in 
Eq. flG.80j) . the absorptive phases cancel, since they do not change sign under charge 
conjugation. The truly CP-violating parts of A however add up, since they change 
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sign for the charge conjugated top decay chain. See the last term of Eq. (1E.67[) . 
which is the CP-violating part of the amplitude squared. 

There is one subtlety in the decay of the neutralino. It decays to equal parts into 
a negatively or positively charged slepton, respectively, leading to the decay chains 

Xi->tf + fe h^x\ + r 2 - (G.84) 
X°->*T+ife %->x?+#. (G.85) 

The two corresponding asymmetries for these two decay modes again differ in sign. 
This is due to a sign change in the spin correlation term of the neutralino T, b D , see 
Eq. (lE.56p . One would have in principle to distinguish the near and the far leptons, 
t\ and £ 2 respectively, not to sum up the asymmetries to zero [U|. But instead an 
asymmetry can be considered of a triple product which includes both leptons from 
the slepton decay, for example, 

T = (pn x Pe 2 ) ■ P^3 = [P/nPfe, Ptsl- (G.86) 

Then a distinction of the near and far leptons £i, £2 from the slepton decay is not 
needed, since there is an additional sign change due to: 

[P/i,P/ a ,P* 3 ] = -[P4>P4,P4s]- ( G - 87 ) 

Thus one only would have to distinguish the charges of the two leptons by considering 
the triple product [64] 

T=\p e +,pe-,p e +], (G.88) 
thereby constructing the CP asymmetry 

A C p = ^A([p e +,Pi~,p e +]) - ^A([pe+,p e -,p £ -]), (G.89) 
which is no longer sensitive to absorptive phases. 
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